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Abstract 

Gravitation has interesting consequences for electromagnetic wave propagation in a 
vacuum. The propagation of plane waves with phase velocity directed opposite to the 
time-averaged Poynting vector is investigated for a generally curved spacetime. Con- 
ditions for such negative-phase-velocity (NPV) propagation are established in terms 
of the spacetime metric components for general and special cases. The negative energy 
density implications of NPV propagation are discussed. 

Keywords: General theory of relativity; negative refraction; negative energy density; plane 
wave propagation. 

PACS: 04.20. Cv, 03.50.De 



1 Introduction 

Negative refraction of a plane wave at the planar interface of two linear, isotropic, homoge- 
neous materials is said to occur when the projections of the real parts of the wave vectors of 
the incident and the refracted plane waves on the interface normal are oppositely directed. 
Then the real part of the wave vector and the time-averaged Poynting vector are parallel 
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in one material, but antiparallel in the other [1]. We call the latter kinds as negative- 
phase-velocity (NPV) materials, but at least two other names have common currency too: 
left-handed materials, and negative-index materials. In order to extend the phenomenon of 
negative refraction to anisotropic materials, NPV materials are characterized by the nega- 
tivity of the projection of the real part of the wave vector on the time-averaged Poynting 
vector. 

Since the beginning of the year 2000 [2], NPV materials have excited much theoreti- 
cal as well as experimental interest. Initial disbelief and criticism in some sections of the 
electromagnetics research community [3, 4] eventually gave way to widespread, but perhaps 
still not universal, acceptance of NPV materials with unequivocal demonstrations by several 
independent groups [5]-[7]. A simplistic expression of the (monochromatic) electromagnetic 
energy density turns out to yield negative values [8], which are generally held as impossible 
in the electromagnetics research community, but more sophisticated investigations indicate 
that the electromagnetic energy density in NPV materials is indeed positive when account 
is taken of the frequency-dependent constitutive properties [9]. 

Perhaps the potentially most useful application of NPV materials is for the so-called 
perfect lenses [10]. Once satisfactorily designed and fabricated, such lenses - although not 
really perfect [11, 12] - could find widespread use in modern optics, for communications, 
entertainment, and data storage and retrieval. More uses would emerge with ongoing re- 
search on anisotropic NPV materials, particularly with negligibly small dissipation in certain 
frequency ranges. 

Instead of concentrating on devices, we turned our attention to the marriage of the special 
and the general theories of relativity (STR and GTR) and NPV propagation of electromag- 
netic fields. We found, a few months ago, that materials that appear to be of the non-NPV 
type to relatively stationary observers can appear to be of the NPV type to observers moving 
with uniform velocity [13]. That result permitted us to envisage STR negative refraction 
being exploited in astronomical scenarios [14] such as for the remote sensing of planetary 
and asteroidal surfaces from space stations. Application to remotely guided, extraterrestrial 
mining and manufacturing industries can also be envisioned. Furthermore, many unusual as- 
tronomical phenomenons would be discovered and/or explained via STR negative refraction 
to interpret data collected via telescopes. 

Ordinary vacuum (i.e., matter-free space) appears the same to all observers moving at 
constant relative velocities. Therefore, NPV propagation in vacuum cannot be observed by 
such observers. This could lead one to believe that NPV propagation is impossible in huge 
expanses of interstellar space. However, gravitational fields from nearby massive objects will 
certainly distort electromagnetic propagation, which is a principal tenet of the GTR and 
is indeed used nowadays in GPS systems, so that NPV propagation under the influence of 
a gravitational field required investigation. In a short communication [15], we derived a 
condition for NPV propagation to occur along a specific direction in a region of spacetime, 
with the assumption of a piecewise uniform but otherwise general spacetime metric. As the 
consequences of such a possibility are highly relevant to further exploration of outer space 
as well as for industrial operations therein, we undertook a more general study, the results 
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of which are being reported here. 

The plan of this paper is as follows: In Section 2, electromagnet ism in generally curved 
spacetime is transformed from a covariant to a noncovariant formalism, wherein vacuum 
resembles a bianisotropic "medium" which enables planewave propagation to be examined 
using standard techniques. A piecewise uniform approximation of the spacetime metric is 
then implemented in Section 3, and thereby a condition for NPV propagation is derived. 
Section 4 is devoted to a discussion of energy density, and the paper concludes with a 
summary in Section 5. 

2 Electromagnetism in Gravitationally Affected Vac- 
uum 

The effect of a gravitational field is captured by the metric g a p which is a function of space- 
time x a and carries the the signature (+, — , — , — ). § In the absence of charges and currents, 
electromagnetic fields obey the covariant Maxwell equations 

fa(3;u + f/3u;a + fua;/3 = , ti* .p = , (1) 

where f a p and h a/3 are, respectively, the covariant and the contravariant electromagnetic 
field tensors whereas the subscript .„ indicates the covariant derivative with respect to the 
vth spacetime coordinate. 

2.1 Noncovariant equations for vacuum 

Following common practice [16]-[18], the Maxwell equations (1) may be expressed in nonco- 
variant form in vacuum as 

fa/3,u + f(3u,a + f va.fi = , 

wherein g = det [g a fs} and the subscript jV denotes ordinary differentiation with respect to 
the isth spacetime coordinate. Although the generalization of the Maxwell equations from 
noncovariant to covariant formulations is not completely unambiguous [19], we adopt the 
standard generalization (1) in the absence of experimental resolution of the ambiguity. 
Introduction of the electromagnetic field vectors 

Ei = feo , 

-Be ( ~L/2,^E£ mn fran 

D e = (-g) 1/2 h e \ 
H e = ( l/2)e tmn {-g) 1/2 H 

§ Greek indexes take the values 0, 1, 2 and 3; Roman indexes take the values 1, 2 and 3; x° — ct where c 
is the speed of light in vacuum in the absence of all gravitational fields; whereas x 1 ' 2 ' 3 are the three spatial 
coordinates. 
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with Ei mn being the three-dimensional Levi-Civita symbol, allows us to state the Maxwell 
equations in the familiar form 



Be,e — , B ifi + e emn E m:n — 
De,e = , — Di$ + et mn H mjn = 



(4) 



The accompanying constitutive relations of vacuum can be written for the electromagnetic 
field vectors as 



B(, r )i m H m ££mn B n 

where 



(5) 



/ \l/2 9 
llm = ~ {-9) 

9oo 

T-, 90m 



(6) 



9oo 

The most important of the foregoing equations can be expressed in SI units as 

V x E(ct,r) + ^- t B(ct,r) = 0, (7) 

VxH(ct,r)-^D(ct,r) = 0, (8) 

Diet, r) = e Met, r) • Eict, r) - - F(ct, r) x Hict, r) , (9) 
= c 

B(ct,r) = Hol(ct,L) 'K(ct,r) + -T_(ct,r) x E(ct,r) , (10) 
= c 

where space r has been separated from t, the scalar constants e and fi denote the per- 
mittivity and permeability of vacuum in the absence of a gravitational field; j(ct,r) is the 

dyadic-equivalent of 7^ m , and T(ct,r) is the vector-equivalent of T m . These four equations 
are stated in the usual style of 3-dimensional vectors and dyadics for convenience, but the 
spacetime is still curved. 

2.2 Partitioning of spacetime 

Let the spacetime region of interest X be partitioned into an appropriate number of subre- 
gions ( m ) X, (m = 1, 2, 3, . . . ). In the mth subregion, the nonuniform metric g a p is written as 
the sum of the uniform metric ^g a p and the nonuniform residual metric ^d a p as follows: 

9af3 = im) 9af3 + {m) d a/3 , (11) 

Notice that, whereas the curved spacetime metric g a p is transformable into the Lorentzian 
metric rj af3 = diag (1,-1,-1,-1) of flat spacetime at every point in X, in accordance with 
the Einstein equivalence principle, the transformation is not universal [20]. That is to say, 
since g a p represents a generally curved spacetime it cannot be replaced by the flat spacetime 
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metric r] af3 everywhere in X. Furthermore, there is no reason for ^g a p to be transformable 
into r\ a p at even one point in ( m ^X. 

The Maxwell curl postulates read as follows in ( m ^X: 



V x E(ct,r) 



d_ 
di 



to M 7 * ^K{ct,r) + - x |W,r) 
= at c at 



to {m U(ct,r) -H(ct,r) + - {m) $(ct,r) x E(ct,r) 



VxH(ct,r) = 



e {m) l ' Qj.E{ct,r) - - c (m) f x -^K{ct,r) 



d 

+ m 



e (m) 0(ct,r) -E(ct,r) - - (m) $(ct,r) x H{ct, : 



(12) 



(13) 



Here, (m) 7 and (m) £ are related to ^g a p, and (m) and (m) $ to (m) d a/3 , in the same way 
that 7 and £ are related to g a g. 



2.3 Piecewise Uniform Approximation 

Equations (12) and (13) are complicated. Therefore, for preliminary analysis, ( m ^g a p can be 
selected appropriately for ^ X and ^ m ^ Q/3 can be ignored. This piecewise uniform approxi- 
mation leads to the simpler equations 



V x E(ct, r) 

V x H(ct, r) 



to (m) l ' §- t H(ct,r) + - c x ^E(ct,r) 

Ho M 7 ' ^E{ct,r) - - x ?-H(ct,r) 
= at c at 



(14) 
(15) 



for electromagnetic fields in ( m ^X. 

The nature of the "medium" implicit in (14) and (15) is worth stating: This medium 
is spatially homogeneous and local, it does not age, and it reacts purely instantaneously - 
just like vacuum in the absence of a gravitational field. However, it is bianisotropic. As 
7 is real symmetric, both the permittivity and the permeability dyadics (i.e., e ^7 and 

fj, ("^7, respectively) are orthorhombic and have the same eigenvectors. Furthermore, the 
gyrotropic-like magnetoelectric terms on the right sides of the two equations can be removed 
in the temporal-frequency domain by a simple transform [21], so that this medium is unire- 
fringent despite its anisotropy. Unless ( m )£ is a null vector, this medium is not reciprocal in 
the Lorentz sense [22]; despite its nonreciprocity in general, the medium satisfies the Post 
constraint [23]. Finally, the medium is nondissipative [24, p. 71]. 



3 Planewave analysis 

In this section, we turn to the propagation of plane waves in a medium characterized by the 
constitutive relations (9) and (10), within the subregion ^ X of X. The presuperscript ^ 
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is not generally used henceforth in this section to avoid cluttering up the equations. It is to 
be understood that all equations in this section hold in with 7 = ("^7 and T_ = ( m )f\ 

We derive an expression for the time-averaged Poynting vector (P) t associated with 
the plane wave with wave vector k, and thereby establish sufficient conditions for NPV 
propagation, as signalled by 

k'{P)t<0. (16) 

3.1 Fourier representation 

The following three-dimensional Fourier representation of the electromagnetic fields is ap- 
propriate for further analysis: 

E(ct,r) = - / / g(uj/c,k) exp[i(k' r- cut)] dujdk 1 dk 2 (17) 

^ J —00 J —00 J — OO 
roc poo poo 

H(ct,r) = - / / H(u/c,k) exp[i(k-r-Lut)] dujdk 1 dk 2 . (18) 

^ J — 00 J — 00 J— 00 

Here, % = y/—T, k = fcif^ + k 2 v 2 + k 3 v 3 is the wave vector with k\ + k\ + k\ = k • k, to is the 
angular frequency, and 

g(to/c,k) = A a (u/c,k)^(u;/c,k) + A b (u/ 0,^^(0; /c,k) , (19) 
H(u/c,k) = A a (uj/c,k) ha(u/c,k) + A b (u / 'c, k) h 6 (u; / c, k) . (20) 

The complex-valued scalars A a b are unknown amplitude functions that can be determined 
from initial and boundary conditions, and the planewave functions e a fe and h a fe are as yet 
unspecified. Clearly, the electromagnetic fields are thus represented in terms of an ensemble 
of propagating plane waves, both propagating (real-valued k 3 ) and evanescent (complex- 
valued /c 3 ). 



3.2 Propagating plane waves 

Further interest being only in propagating waves, we set k 3 G R. The determination of 
e a b and h a b follows the same path as for propagation in a simply moving medium that is 
isotropic dielectric-magnetic at rest [24, chap 8]. 

Substituting (17) and (18) in (14) and (15), we find that 

pxe a b (u/c,k) = ujn g • h a b (a;/c, k) , (21) 

= -^ e o2 •§ a)6 (u;/c,fc) , (22) 

(23) 

where 

p = k--T. (24) 
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Hence, 



while 



ha,6(«Vc,£) = —7 1 * \p x e ajb (a;/c, fc)] , 
|(pxl)- (adj£) '(px|) + Q J III} '^,6 = 0, 



(25) 



(26) 



where "adj" stands for the adjoint, J is the identity dyadic, and \^\ is the determinant of 7. 
As 7 is symmetric, the foregoing equation can be further simplified to 



(~) \l\-P'l'P 

Equation (27) has nontrivial solutions if 

(~) \l\-P'l'P 



L + vp ' 7 



l + pp '7 



= 0. 



(27) 



(2? 



From this condition, the dispersion relation 

P'l'P- (^) 







(29) 



to determine ks for specific {u, k±, A^} emerges. In the three-dimensional p-space, this 
relation represents the surface of an ellipsoid [25, Sec. 3.5.4]. 
Substituting (29) in (27), we obtain 

PP*l*e a)b = 0, (30) 

whence 

P ' J ' e a , 6 = . (31) 

Thus, both and e b must be orthogonal to p • 7, but neither of the two is generally orthog- 
onal to p. A similar exercise yields 

P'l'K, b = 0, (32) 

so that both h a and h fe must be orthogonal to p • 7 but not necessarily to p. 

The selection of e a fe to satisfy (31) produces a trilemma, which can be explained as 
follows: We can always choose two unit vectors w and y that are orthogonal to each other 
as well as to p. Then, without loss of generality, (31) is satisfied by 



7" 



w 



r 1 



(w + qy) 



*a = fr=i ; .,.1 , & 



1 7 * Ul\ 



ir 1 



(w + qy)\ '' 



(33) 



where gel, while h afe can be obtained from (25). The following three conditions appear 
reasonable in order to fix q: 
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(i) § a • e b = 0, 

(ii) ha ' hj, — 0, and 

(iii) e a x h b = (or, equivalently, e 6 x = 0). 

In general, the three conditions turn out to be mutually exclusive, i.e., only one of the three 
can be fulfilled. We chose 

w + qy = pxe a (34) 
in order to fulfil the third condition; thus, 

g 1 -W J" 1 * (P X £») 

5a = TT— i | , e 6 = — — f -; r— . (35) 

3.3 Negative— phase— velocity propagation 

The rate of energy flow for a specific {co,k} in ^ X is obtained by averaging the Poynting 
vector over one cycle in time; thus, using (25) and (35) we find that 

(P>t = 7T^fl A «| 2 5a-7-e a + |A 6 | 2 e 6 .7-eJ7-p- (36) 
Zunoljj 1 = = j = - 

The important result therefore emerges: the NPV condition (16) is satisfied provided that 

[\A a \ 2 e a • J • e a + \A b \ 2 e b • g • e 6 ] k • J • p < . (37) 

Let us elaborate on the NPV inequality (37) for positive- and negative-definite 7 and 
indefinite 7 [27]. 

3.3.1 Positive- and negative-definite 7 

Since 7 is a 3x3 dyadic, we have ^ • 7 • e a j /|7| > and (^e b • 7 • e 6 ^ /|7| > when 7 

is either positive-definite or negative-definite. Thus, the NPV inequality (37) immediately 
reduces to 

k-^'p<0. (38) 

We therefore observe that NPV propagation arises for negative-definite 7 if the gyrotropic- 

like magnetoelectric term T_ = 0, whereas NPV propagation is not possible for positive- 
definite 7 if r = 0. 
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In order to establish the sign of k • (P) t , let us introduce the mutually orthogonal basis 
vectors b x ,b_ 2 and 63, where b 3 is parallel to 7 • £ but b x and b 2 lie in the plane perpendicular 

to 7 • f. With respect to these basis vectors, we express k and f as 



k = Kib x + /t 2 6 2 + ^363 
f = Giftj + G 2 b 2 + G 3 b 3 



(39) 



The definiteness of 7 ensures that f is not perpendicular to 7 • f . Thus, G 3 7^ and we 
have 

* = ft + (*i ~ S«s j &! + ^2 - Sk 3 ) k 2 . (40) 

Equivalently, the wave vector k may be written in the form 

k = Pl -T + p 2 z, (41) 
c 

with pi j2 G 1R being scalar constants and z being a unit vector in the plane perpendicular to 
7 • f. It follows that 

k'J_'P~{~) P 1 ~ £ * T 'tl + plz *7 (42) 

Hence, NPV propagation is a consequence of the inequality 

/to \ ^ ~ ~ 

pk*i*^< (-) Pi(i-Pi)r-j-r (43) 

being satisfied. 

We emphasize that the NPV condition (43) applies for an arbitrarily oriented wave vector 
k. Two particular cases are worthy of special mention. First, if k lies in the plane perpen- 
dicular to 7 • £ (i.e., p\ = 0), then NPV propagation cannot occur regardless of the value of 

p 2 or orientation of z. Second, suppose the wave vector k is aligned with T_ (i.e., p 2 — 0). 
Then the NPV inequality (43) is satisfied for all Pl G (0, 1). 

3.3.2 Indefinite 7 

The dyadic 7 need not be positive- or negative-definite. For example, the dyadic 7 corre- 
sponding to the Kerr metric of a rotating black hole is indefinite in the ergosphere region 
[20]. Where 7 is indefinite, the sufficient conditions 



for NPV propagation emerge from (37). 



^(e 6 *2 '§b) (k *1 *p) < 



(44) 
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3.4 Summary of main results 



Let us summarize our main results in this section: For the subregion of generally curved 
spacetime (™)X E X, with curvature specified by the uniform metric ^g a /3, plane waves 
propagate with phase velocity directed opposite to the direction of the time-averaged Poynt- 
ing vector provided that the inequality (37) holds. When ^7 is either positive- or negative- 
definite, the NPV inequality (37) reduces to the simpler relation (43), whereas the sufficient 
conditions (44) apply when ^7 is indefinite. 



4 Energy density 

When dealing with plane waves in linear, homogeneous materials, it is common to define the 
time-averaged electric and magnetic energy densities as 



(\N e (u/c,k,r)) t = \Re[E(u/c,k) • D*(u/c, k)} cxp [-2Im(fc • r)] 
(\N m (u/c,h,r))t = |Re[H>/c,A) • B(w/c, k)] exp [-2Im(fc • r)] 



(45) 



where the asterisk indicates the complex conjugate, while D(uj / c, fc) and B(u/c, k) are defined 
similarly to E(u>/c,k) in (17). According to the Maxwell curl equations, 



k x E(u/c,k) = uj B(u/c, fc) 
k x H(u/c,k) = -cuD(cu/c,k) 

therefore, 

(\N e (uj/c,k,r)) t = ±Re{k* • [E(u/c,k) x H>/c, A;)]}exp [-2 Im (fc • r)] 
(\N m (u;/c,k,r)) t = ±Re{k • [E(w/c,fc) x H>/c, £)]} exp [-2 Im (fc • r)] 

The total time-averaged electromagnetic energy density (W(ou/c,k,r)) t is the sum 



(46) 



(47) 



(W(oj/c,k,r)) t = (\N e (u/c,k,r)) t + (\N m (u;/c,k,r)) t 

= Re (fc) • Re ]E(u/c, fc) x H* (w/c, fc)] exp [-2 Im (fc • r)] . (48) 

As the time-averaged Poynting vector 

(P(cu/c,k,r)) t = - [E(u/c,k) x H*(cu/c,k)} exp [-2 Im (fc • r)] , (49) 

it follows that 

(W(w/c, fc, r))t = - Re (fc) • (P(w/c, fc, r)) t exp [-2 Im (k • r)] . (50) 

Thus, the electromagnetic energy density (as calculated in this paragraph) associated with 
a NPV plane wave must be negative. 
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Relevant to Section 3, 

D(u/c,k) = e„7 •E(iv/c,k) - -f x H(w/c,k) 
= c 

B(u/c,k) — /i 7 * H(u;/c, fc) + - f x E(u/c,k) 
= c 

and the possibility of negative (\N(cu/c, k, r)) t for propagating plane waves in gravitationally 
affected vacuum emerges. 

The possibility of a negative electromagnetic energy density requires discussion. In the 
research on isotropic, homogeneous, dielectric-magnetic NPV materials, the negative value 
has been noted [8]. Equally important is the fact that such materials have been artificially 
fabricated as composite materials comprising various types of electrically small inclusions, 
and their planewave response characteristics (over limited u;-ranges) are substantially as 
predicted [7]. This means implies the aforementioned procedure to compute (W(o>/c, k, r))t 
may not be always correct. Indeed it is not, because it applies only to nondissipative and 
nondispersive mediums. When account is taken of the dissipative and the dispersive nature 
of the NPV materials [26], (\N(u/c, k, r)) t does turn out to be positive [9]. 

However, the medium in Section 3 is nondissipative and nondispersive, so that the fore- 
going paragraph does not apply — but it does provide the basis for the following argument. 
Electromagnetic energy densities for plane waves, howsoever computed, are not necessarily 
indicative of the true picture. This is because an electromagnetic signal is of finite spatiotem- 
poral extent, while plane waves are infinitely extended over the entire spacetime; indeed, it 
can be argued that a plane wave has infinite energy! Therefore, the energy density of a signal 
is meaningful, but the time-averaged energy density of a plane wave may not be. In com- 
puting the energy density of a signal, one must consider the bandwidth in the uo © k domain. 
Since the NPV conditions in Section 3 appear unaffected by u but not by the direction of 
propagation, NPV plane waves could appear in gravitationally affected vacuum as part of 
a pulsed electromagnetic beam (of finite cross-section) which has positive and finite energy 
density. 

A proposal to overcome the negative value of {W(u/c,k,r))t in NPV materials is to 
fabricate them out of active inclusions [28]. Whereas passive inclusions are modeled in 
terms of resistances, capacitances, and inductances, the modeling of active inclusions invokes 
amplifiers as well. In other words, there is a source of energy to offset negative (\N(u/c,k, r)) t . 

Reverting to (11), we see that ^d a /3 could be considered as a spatiotemporally non- 
homogeneous source term. The effect of this term must be included in all energy density 
calculations, in addition to the effect of the finite spatiotemporal extent of any electromag- 
netic signal. In other words, the spatiotemporal fluctuations of gravitation can act as a 
source term. Thus, one must consider the total energy density, not just the electromagnetic 
part of it. 

Astrophysics researchers have formulated several different energy conditions for classical 
(i.e., nonquantum) GTR: all are just conjectures lacking rigorous proofs from fundamental 
principles and were set up simply to prove certain theorems [29]. Violations of these energy 
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conditions are known [30, 31], and negative energy densities are invoked for the formation of 
certain black holes [32] as well as for the phenomenon of superradiant scattering of electro- 
magnetic, gravitational, and scalar waves [33, Sec. 12.4]. In fact, two astrophysicists have 
recently written [34]: 

It is often (mistakenly) believed that every kind of matter, on scales for which 
we do not need to consider its quantum features, has an energy density that is 
everywhere positive. 

This situation arises because the local energy density of a gravitational field cannot be 
defined uniquely in GTR. However, the notion of the total gravitational energy-momentum 
of an isolated system — such as ADM energy-momentum [35] - - is available in an asymp- 
totically flat spacetime; see also [36, 37]. But there is no guarantee that ADM total energy 
should be positive. The condition of positivity of energy can only be expected to hold if the 
spacetime is nonsingular and this condition is imposed on matter distribution [38] , and that 
for isolated systems [39, 40]. 

Under certain circumstances, many exotic solutions of general relativity have been shown 
to have negative energy densities. Such studies have exploited the use of quantum fields as 
possible sources of negative energy densities [41, 42, 43]. Unlike classical physics, quantum 
physics does not restrict energy density to have unboundedly negative values (though there 
are some bounds that constrain their duration and magnitude [41, 42]), which then enable 
the quantum fields to be used to produce macroscopic effects. 

In summary, the issue of energy density remains to be carefully investigated for electro- 
magnetic fields in gravitationally affected vacuum, regardless of the satisfaction of the NPV 
condition (43). This will require numerical studies with specific spacetime metrics. A similar 
resolution is needed for the Casimir effect [44, 45]. 

5 Concluding remarks 

We have investigated the propagation of electromagnetic plane waves in a generally curved 
spacetime. Sufficient conditions (37), (43) and (44) for negative-phase-velocity propagation 
are established in terms of the spacetime metric components. We conjecture that research on 
new phenomenons encountered during space exploration — e.g., the anomalous acceleration 
of Pioneer 10 currently being observed [46] - may benefit from NPV considerations. The 
negative energy density implications of NPV propagation require further investigation. 

Acknowledgement: SS acknowledges EPSRC for support under grant GR/S60631/01. 
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